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ABSTRACT 

We present long-baseline interferometry of the Kepler exoplanet host star HD 179070 
(Kepler-21) using the PAVO beam combiner at the CHARA Array. The visibility 
data are consistent with a single star and exclude stellar companions at separations 
~l-1000mas (~ 0.1-113 AU) and contrasts < 3.5 ma gnitudes. Thi s resul t supports 
the validation of the 1.6 exoplanet Kepler-21b by iHowell et al.l (|2012f) and com- 
plements the constraints set by adaptive optics imaging, speckle interferometry, and 
radial velocity observations to rule out false-positives due to stellar companions. We 
conclude that long-baseline interferometry has strong potential to validate transiting 
extrasolar planets, particularly for future projects aimed at brighter stars and for host 
stars where radial velocity follow-up is not available. 

Key words: stars: individual: HD 179070 - planets and satellites: individual: Kepler- 
21b - techniques: interferometric. 



1 INTRODUCTION 

The NASA Kepler Mission aims to find extrasolar planets 
in the habitable zones of solar-type stars through the detec- 
tion of brightness dips as planets cross the stellar disc. While 
Kepler has been highly successful in finding exoplanet can- 
didates, ground-based follow-up observations are important 
to confirm the detections. The most common astrophysical 
false positives for Kepler involve stellar companions that re- 
main unresolved due to Kepler's large pixel size (~ 4"). False 
positives can be divided into companions that are physically 
bound to the target star (hierarchical triple systems) and 
companions that are either in the foreground or the back- 
ground of the target due to chance alignment (blends). In 
both transit-like shape can be mimicked by eclipses 
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of a stellar companion or transits of a planet around the 
secondary companion. 

For large (Jupiter- and Neptune-sized) planets, candi- 
dates can often be confirmed using radial velocity obser- 
vatio ns, giving a direct e s timate of the planet's mass (see, 
e.g. iBorucki et all l20ld; I Koch et al. | l2010l ; lLatham et al,l 



l20ld : Cochran et al.l 201 if ), while transit timing variations 



can be used to confirm p l anets in multiple systems (see, 
e.g., iFabrvckv et al. 1 12012| ; iFord et all 120121 : ISteffen et ail 
2012). For many super-Earths and Earth-sized planets, how- 
ever, the Doppler signature is typically too small com- 
pared to the intrinsic stellar variability, and transit tim- 
ing variations might not be detected. In these cases, can- 
didates are validated by excluding as many false-positive 
scenarios as possible. The first stage in this process uses 
the Kepler data to detect signatures of stellar compan- 
ions through photocenter shifts and the comparison of 
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trans it depths (see, e.g.. iBatalha et all 120101 : iBrvson et alJ 
|2010| ). followed by st atistical modeling o f potential blending 
scenarios (see, e.g. , Torres et al.l l201ll ; iFressin et al.l l201ll ; 
iMorton fc Johnson! l201lF 7 These constraints are then com- 
bined with ground-based follow-up observations such as 
spectroscopy, spec kle interferometry and adaptive optic s 
imaging (see, e.g.. iGautier et al]|2010l ; iHowell et al.ll201ll) . 
High-angular-resolution observations using long-baseline in- 
terferometry offer a powerful tool to complement these 
methods and extend the parameter range that can be ex- 
cluded, particularly for close-in (both bound and unbound) 
com panions. 

(Howell et all ||2012T ) recently reported the detection of 
Kepler-21b, a 1.6 R(q planet in a 2.8 d orbit around the 
bright (V = 8.3) F6IV star HD 179070. Extensive follow- 
up observations have been used to rule out a false positive 
detection in this system. Speckle interferometry rules out 
any unbound or bound stellar companions at separations 
> 0.2" (> 22.6 AU) for contrasts up to 5 mag, and at sepa- 
rations > 0.05" (> 5.6 AU) for contrasts up to 4 mag. Adap- 
tive optics imaging revealed a ~14.5mag companion at a 
separation of 0.7" which, however, was found unable to mim- 
ick the observed transit shape due to its low mass. Radial- 
velocity time series were also obtained, and showed no sig- 
nificant variation over 85 days at a 5.6 m/s level. Putting 
all these constraints together, the only possibilities escap- 
ing direct detection are close-in bound companions within 
~4mag of HD 179070 in nearly face-on orbits (causing no 
detectable RV signature) with a similar radial velocity to 
the target star, and close-in blends within ~ 7 mag that may 
have escaped spectroscopic detection. Here, we present long- 
baseline interferometry obser vations to compleme nt and ex- 
tend the validation efforts bv lHowell et al.l (|20121 i. 



Table 1. Spectral types, photometry and expected angular sizes 
of stars in this study. Brackets indicate the last two digits of the 
1<t uncertainties. 



HD 


Sp.T. 


V 


K 


9v-K (mas) 


179070 


F6IV 


8.262(11) 


6.945(18) 


0.169(09) 


174260 

178591 1 

183204 


B8V 
B5V 
AOV 


7.323(05) 
7.130(57) 
7.425(07) 


7.465(18) 
7.191(24) 
7.386(21) 


0.103(05) 
0.119(06) 
0.110(06) 



1 HD 178591 is an ellipsoidal variable with an amplitude of 
~ 0.05 mag. 
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2 OBSERVATIONS 

We have observed HD 179070 as part of our interferometric 
follow-up campaign of Kepler stars using the PAVO (Pre- 
cision Astronomica l Visible Observations) beam combiner 
(|lreland et al.l 120081 ) at the CHARA (Center for High An- 
gular Resolution Astronomy) Array (|ten Brummelaar et al.l 
120051 ). PAVO is a three-beam pupil-plane beam combiner 
optimized for high sensitivity and angular resolution, record- 
ing visibilities over a spectral bandpass of ~ 650-800 nm 
with a limiting magnitude in typical seeing conditions 
of R < 8 mag. Using baselines reaching up to 330 m, 
PAVO/CHARA is capable of resolving angular sizes down 
to ~ 0.3 mas. For mo re details on the inst rument and data 
reduction, we refer to llreland et all l|2008h . 

HD 179070 was observed on 2 July 2011 in two-telescope 
mode using the Sl-Wl (278 m) baseline in excellent see- 
ing conditions. Two scans were obtained, which were inter- 
leaved with observations of three different calibrator stars. 
Using various catalogs available in the literature, calibra- 
tors are typically chosen to be single field stars in close 
vicinity (< 10°) to the target star. Table 1 lists the spec- 
tral types, photometric properties, and expected sizes of all 
stars in this study. The predicted si zes were calculated us - 
ing the (V — K) relation given by iKervella et all l|2004l ). 
V magnitudes have been extracted from the Tycho catalog 
l|Perrvman fe E SA 1997) and transformed into the Johnson 
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Figure 1. Raw visibility versus wavelength for the target 
HD 179070 (top panel) and the calibrator stars (bottom panel). 



system using the calibration of lBesselll {2000). K magnitudes 
have been obtained from the 2MASS catalog (|Cutri et al.l 
l2003l;ISkrutskie et al.ll2006t ). We have tested the photometry 
for reddening by comparing the observed (B — V) colors with 
a list of intrinsi c color s as a function of spectral type given by 
ISchmidt- Kaleij(|l982l ). The observed colors have been found 
to be compatible with the spectral types for all stars except 
for the calibrator HD 178591, which is classified as an ellip- 
soidal variable by Hipparcos. We have accounted for this by 
adding a systematic uncertainty corresponding to the ob- 
served variability amplitude (~ 0.05 mag) to the statistical 
photometric uncertainty in the V band. Note that the po- 
tential companion of HD 178591 can be expected to have 
a negligible influence on the interferometric measurements. 
The final uncertainties for the predicted angular diameters 
were calculated by adding a conservative 5% calibration un- 
certainty in quadrature to the photometric uncertainty. 



© 0000 RAS, MNRAS 000, 000-000 



Validation of Kepler-21b using PAVO/CHARA 3 



1.4 




Figure 2. Calibrated squared visibility versus spatial frequency for HD 179070. The red solid line shows the best fitting single disc model 
with a diameter of $ld = 0.13 ± 0.02 mas. The areas marked by blue dashed, dashed-dotted and dashed-triple-dotted lines show the 
range of minimum squared visibilities expected for companions > 1 mas with contrasts of 1.5, 3 and 4.5 magnitudes, respectively. The 
inset shows a close-up of the data. 



3 RESULTS AND DISCUSSION 

Figure 1 shows the raw squared visibility measurements 
across the PAVO passband for the target and the calibra- 
tors. As expected from their predicted sizes, the visibility 
levels of the calibrators are very similar. Note that the raw 
measurements are considerably lower than the expected vis- 
ibilities (~ 0.9) due to atmospheric turbulence and opti- 
cal aberrations. HD 179070 remains practically unresolved 
in our observations, with raw visibilities at very similar lev- 
els to the calibrator stars, and no significant visibility change 
as a function of wavelength. More importantly, the similarity 
between target and calibrator scans shows that the visibil- 
ity curve of HD 179070 is consistent with a single star. For 
any companion for which the interference patterns (fringe 
packets) of the two stars overlap in delay space, a periodic 
visibility modulation would be observed, while any incoher- 
ent flux from companions at larger separatio ns would cause 
a drop in the observe d visibility (see, e.g.. iMonnierl I20fj3 : 
ten Brummelaanl2007l '). The interferometric field of view for 
each case depends on the coheren ce length and the projected 
baseline ()ten Brummelaarlll995l ). and for our observations 
corresponds to <30mas and 30-1000 mas, respectively. 

To illustrate this more clearly, Figure 2 shows the cal- 
ibrated squared visibility data of HD 179070 as a function 
of spatial frequency. Visibilities were calibrated by dividing 
the calibrator data by the predicted sizes to obtain a system 
visibility, which was then used to correct the target visibil- 
ity. We then fitted th e limb -darkened disc model given by 
lHanburv Brown et al.l (1 19741 ) to th e data using the method 



2011) to the stellar parameters presented bv lHowell et al.l 



2012). The red line in Figure 2 shows the best-fitting single- 



described by Derekas et al. (|201ll ). We used a linear limb- 
darkening coefficient for the i?-band of = 0.5197, derived 
from the closest matching grid point of IClaret fc Bloemenl 



disc model, yielding a diameter of (9ld = 0.13±0.02 mas with 
a reduced x 2 — 1.4. This diameter agrees with the diameter 
of 0.15 ± 0.01 mas co nstrained from the asteroseismic radius 
(R/Rq = 1.86±0.02. lHowell et al.l2012T) and Hi pparcos par- 
allax (n = 8.86 ± 0.58 mas, van Leeuwenl [20071 ) . The areas 
marked by blue dashed, dashed-dotted, and dashed-triple- 
dotted lines illustrate the minimum squared visibilities ex- 
pected for a stellar companion with a contrast of 1.5, 3 and 
4.5 magnitudes compared to HD 179070. Deep minima cor- 
respond to close-in companions (1-30 mas) which will show 
a periodic variation across the PAVO passband, while shal- 
lower minima correspond to wide companions (> 30 mas), 
for which the variation is unresolved within the spectral reso- 
lution of PAVO and hence an overall drop in visibility would 
be observed. Figure 2 suggests that the PAVO observations 
rule out any stellar companions at contrasts of < 3 magni- 
tudes. Note that for stellar companions at even closer sep- 
arations (< 0.1 AU), the PAVO data would solely exclude 
companions along the baseline vector at the time of obser- 
vation, which in this case spans only one epoch. We also note 
that the faint ~ 14 .5 mag companion at ~0.7" detected by 
I Howell et al.l (|2012l ') has negligible influence on our measure- 
ments. 

To establish the magnitude limit more quantitatively, 
we performed 10 s simulations as follows. For each iteration, 
we used the spatial frequencies of our observations to gen- 
erate a synthetic binary model consisting of a primary with 
the expected size of HD 179070 and an unresolved secondary, 
with a separation and contrast drawn from uniform distribu- 
tions between 1-50 mas and 1-6 mag, respectively. We then 
added white noise to each data point with a standard devi- 
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ation corresponding to our estimated relative measurement 
uncertainties. For each simulated dataset, we then fitted a 
binary model to the data and compared the \ 2 value to the 
one calculated from a single disc model with the expected 
size of HD 179070. The simulations showed that in > 99% 
of all cases, the binary model yielded a significantly (> 3 a) 
better fit for contrasts below 3.6 magnitudes. This limit has 
been found to be only weakly dependent on the separation, 
dropping to ~3.4mag for wide (50-1000 mas) binaries. To 
confirm the limit, we performed a Bayesian model compar- 
ison by calculating the odds ratio between both models. 
Equal prior probability of each model was assumed (hence 
simplifying the problem to the calculation of the Bayes fac- 
tor) . The marginal likelihood (evide nce) for each m odel was 
calculated using Nested Sampling l|Skillindl2003 ). We as- 
sumed a Gaussian prior probability for the primary diame- 
ter, and a uniform prior for the separation and contrast of 
the binary model. Due to computational reasons the calcu- 
lation was performed for 150 simulations and restricted to a 
smaller range of close-in (1-10 mas) companions. The com- 
puted Bayes factors consistently favored the binary model 
over the single star model for contrasts < 3.5 mag. As ex- 
pected, this value is more conservative than the limit in- 
ferred from the likelihood ratio, since the Bayesian model 
comparison penalizes the binary model for its added com- 
plexity. Based on these results, we conclude that our data 
would have revealed a stellar companion at separations ~ 1— 
1000 mas (~ 0.1-113 AU) and contrasts < 3.5 magnitudes. 

We note that the above simulations assume that our 
measurement uncertainties are well characterized. Previous 
PAVO observations of the same star over multiple nights 
showed that night-to-night variations in V 2 are at the 2-3% 
level (Huber et al. 2012, in preparation), while the ~ 5 % un- 
certainties in the calibra tor sizes (see Table 1) trans late to 
a 1 % uncertainy in V 2 (|van Belle fc van Belli 120051 "). Both 
these contributions are considerably smaller than our mea- 
surement uncertainties (estimated from the scatter of indi- 
vidual data frames integrated over each scan), which are on 
average 13%. Combined with the low reduced \ 2 of our fit, 
we therefore argue that our measurement uncertainties are 
a good estimate of the total uncertainty in our data. 

The results p resented here c omple ment and extend the 
constraints set by iHowell et ail {2012). For the case of un- 
bound companions, the PAVO limit of < 3.5 mag covers only 
a small range of the possible false-posi tive scenarios, which 
extend down to < 7 mag. Additionally, IHowell et al.l (|2012h 
demonstrated that the probability of a chance alignment of 
a star able to reproduce the transit shape at such close sep- 
arations is very small. Hence, the added information from 
PAVO for unbound companions is negligible. For the case 
of bound companions, Figure 3 illustrates the constraints 
from PAVO and speckle interferometry on the mass and 
separation of a possible secondary. Magnitude limits were 
converted into mass li mits by interpolating a 3 Gyr solar- 
metallicity isochrone bv lBaraffe et al.l fl998), which roughly 
matches the age of HD 179070 as determ ined by the aster- 
oseismic analysis in IHowell et al.l (|2012l ). Additionally, we 
plot the 3-cr detection limit from radial velocity follow-up 
for different inclinations, assuming a circular orbit for the 
hypothetical companion and a simple linear velocity change 
with time. The diagram shows that PAVO rules out face-on 
orbits for bound companions at separations < 5.6 AU, which 
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Figure 3. Constraints on the mass and separation of a sec- 
ondary companion to HD 179070 from PAVO (red dashed line) 
and Speckle interferometry (blue dashed-triple-dottcd line). Mag- 
nitude limits have been converted into secondary masses by in- 
terpol ating a 3 Gyr solar-metallicity isochrone by iBaraffe et alj 
(1998). Black solid lines show the limit from RV follow-up for 
different inclinations, with any companion to the left of the lines 
being detectable with an average radial velocity signal > 16.8 m/s 
(> 3<r) over a timespan of 85 days. Note that any false-positive 
detection for secondaries with masses < 0.8 Mq (marked by a 
long-dashed li ne) can be ruled ou t from the observed Kepler tran- 
sit shape (see IHowell et al]|2012t ). 



could previously not be excluded, and confirms the con- 
straints set by complementary follow-up methods. Since a 
planet transiting a secondary with a mass lower than 0.8 Mr? 
can b e excluded from the Kepler transit shape l|Howell et al.l 
120121 ), the combined constraints virtually rule out any possi- 
ble false-positive scenario involving a gravitationally bound 
companion around HD 179070. 



4 CONCLUSIONS 

We have presented high-angular-resolution observations of 
the exoplanet host star HD 179070 (Kepler-21) using the 
PAVO beam combiner at the CHARA Array. Our data 
clearly rule out stellar companions at separations between 
~ 1-1000 mas (~ 0.1-113 AU) with contrasts of < 3.5 magni- 
tu des. This complemen ts and extends the validation efforts 
by IHowell et al.l (|2012l ). and supports the conclusion that 
the detected transit is due to a 1.6 extrasolar planet in 
an orbit around HD 179070. 

The results shown here demonstrate the potential 
of PAVO/CHARA to validate transiting exoplanet can- 
didates, and complement the existing efforts using long- 
baseline i nterferometry to characterize exoplanet host stars 



(s 



Baines et al] 120091 ; Ivan Belle fc von Braunl 12009; 



Ivon Braun et aL 2011). Using a recent compilation of de- 
tected exoplanets in the NASA Exoplanet Archiv^H we esti- 
mate about half a dozen host stars with transiting exoplan- 

1 http:/ /exoplanetarchive. ipac.caltech.edu/index. html 
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ets to be accessible to observations with PAVO/CHARA. 
Furthermore, there will be a considerable overlap with the 
target sample o f the planned Tran siting Exoplanet Survey 
Satellite fTESS. IPdcker et aljlioosl ). which is aimed at find- 
ing planets around nearby (V < 12) stars. While the contri- 
bution of PAVO to the validation effort of Kepler-2f b is rel- 
atively modest, it can be expected that long-baseline inter- 
ferometry will play a significant role in validating transiting 
extrasolar planets, in particular for future missions aimed 
at bright stars and for cases where precise radial-velocity 
follow-up may not be available. 
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